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In this paper we investigate weak decays of heavy mesons in the framework of a covariant quark 
model, which is based on the Bethe-Salpeter equation in instantaneous approximation. Apart from 
a phenomenological confinement potential, a residual interaction induced by instantons is adopted. 
Masses and many decay observables of light mesons have already been described successfully in this 
model . An appropiate extension allows a unified description of light and heavy systems. 
. Using a set of parameters which are fixed by the mass spectra, we evaluate the form factors of 

semileptonic decays of charmed and bottom mesons. In the heavy quark limit these can be reduced 
5_l to the Isgur-Wise function, which is calculated. Finally the form factors are used to determine the 

non-lcptonic decay rates of B mesons in the factorization approximation. 
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I. INTRODUCTION 



\ In the last few years new and improved data on the spectra and the decays of charmed and bottom mesons have 
s * . become available. The observations of the D\ and the radial excited D*' and B' have been published recently (see 
p3| and references therein). The dominant semileptonic decays of B, D and D s mesons have been measured with 
good precision in the meantime, and data for double Cabibbo-suppressed channels are also available. Yet many new 
results will be provided by the B-factories BaBar, Belle, Hera-B and LHC-B within the next years. 



For the theoretical description of these masses and decays various ansatze are pursued. Lattice gauge theory gives 
good results for the transition form factors for high momentum transfer q 2 , while QCD sum rules are suitable for 
the low q 2 regime. Heavy quark effective theory (HQET), which is based on a systematic expansion in the inverse 
f> |j mass of the heavy quark, provides some model independent predictions, i.e. approximate degenerated mass doublets 
^ . according to different orientation of the heavy quark spin, and a connection between the semileptonic partial decay 
' rate of the pseudoscalar and pseudovector decay channels, which has been experimentally confirmed. Unfortunately 
r"| , it can not predict the masses and decay rates itself, and corrections due to the finite quark mass, at least for the 
L| ■ charm quark, are expected to be substantial. Hence for a consistent description of meson masses and decays over the 
. ,-h \ full kinematic region, constituent quark models, even if the connection to the underlying theory is not quite clear, are 

■ still the most successful tool. 

■ In previous papers we have developed a rclativistic constituent quark model for light mesons, which is based on the 
Bethe-Salpeter equation in instantaneous approximation. Apart from a phenomenological confinement potential we 
adopt a residual interaction induced by instantons. In this model, a very good description of the light meson masses, 
from the ground state nonet up to highest angular momenta, has been achieved. Also many decay observables have 
been calculated in reasonable agreement with the experimental data (sec Q| for a recent update). Motivated by this 
success, the model has been extended for heavy flavours ||. Thereto we do not apply the One Gluon Exchange, 
which is known to give a good description of heavy quarkonia and even of the whole meson spectrum with a suitable 
inclusion of relativistic effects JlT[ ] . Instead we generalize the instanton induced interaction in a naive way^. This 
is done in order to keep the model simple and unified: the parametrization of confinement should be valid for all 
flavours, and only one residual interaction is introduced. Also the question should be raised, when and how the model 
for light mesons fails, if higher quark and meson masses are involved. It turns out, that heavy-light mesons and, to a 
certain extent, heavy quarkonia can be described in that way. The resulting spectra together with a brief review of our 
model are discussed in section II. Knowing the meson amplitudes we then calculate the semileptonic decays without 
further parameters in section III. This is done in order to test these amplitudes rather than to determine decay rates 



: For a similar investigation of weak meson decays in the Bethe Salpeter framework with the One Gluon Exchange see Q. 
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or CKM matrix elements with high precision. Therefore we concentrate on the dominant processes B -> D^tv and 
on the CKM-favoured decay channels of D and D s mesons. Finally in section IV we investigate non-leptonic decays 
of B mesons, and section V contains our summary. 



II. THE BETHE-SALPETER MODEL 



Since the model has been described earlier in greater detail |^||, we will only briefly review the main features. The 
model is based on the Bcthc-Salpeter equation for qq bound states 



x p (p) 



V) 



(1) 



for the Bethe-Salpeter amplitude x P (p) := (0\Tip(p)'4>(p)\P) ■ Here \P) denotes a bound state with mass M and total 
momentum P, P 2 = AI 2 . In our ansatz the full quark propagators Sf are approximated by free fermion propagators 
Sf(p) « ij (rf — m,i + ie), where the constituent quark masses m, are introduced, which are treated as free parameters 
of the model. Furthermore, the irreducible interaction kernel K is assumed to be instantaneous in the rest frame of 
the meson, K(P;p,p')\ p _, M gs = V(p,p'). These assumptions lead to the (full) Salpeter equation 



$(p) = Ar(#7° 



dV V{p,p')$(p') 



(2tt) 3 M + UJI+ LU 2 



7°A+(-p)-A+(p) 7 



dV V(p,p')1>(p') 



(2tt) 3 M-uJi- uj 2 



7°A 2 -(-p) 



(2) 



for the Salpeter amplitude $(p) := J ^~X P (p)\p = im o) ■ Here u>i := \J p 2 + mf and A* denotes the projector on 
positive and negative energy solutions of the Dirac equation. 

As ansatz for the interaction kernel we apply the sum of a phenomcnologically motivated confinement potential and 
a residual interaction, which is induced by instantons. Confinement is parametrized by a linearly rising potential in 
configuration space with an adequate Dirac structure, symbolically written as 

V C (r) = (a c + b c ■ r) T <E) T. 

To estimate the influence of the Dirac structure, two possibilities are taken into account, in the following referred 
to as model A and model B: a combination of scalar and time-like vector, which has already been discussed in [|]||] 
and which is known to minimize spin-orbit splittings, and a chirally invariant combination of scalar, pseudoscalar and 
vector type, previously investigated by (8) and || . 

The additional instanton induced residual interaction is based on the work of 't Hooft [[[o| and Shiftman et al. [ pr| . 
The effect of QCD instantons on the quarks leads to an effective two body interaction, which can be expressed by the 
following kernel in momentum space 

j-^Vu^p'Mp') = 4G J^R A (p,p<) • (tr 1 + tr [ 7 5 <J>(p")] 7 5 ) ■ (3) 



The strengths of the coupling are collected in the flavour coupling matrix 



-3/1/2 : A = h + h = h 

G hh,Sifi = { 3/1/3 : A = h ^ h = h 
: otherwise, 



(4) 



with flavour indices The originally point like interaction is regularized by the function i?A, for which a gaussian 
form is used. 

Heavy mesons have been included in this model to achieve a unified description of all mesons. Especially the 
parametrization of confinement should be universal. As shown in ||, this is done most successfully by extending 
the instanton induced residual interaction in a naive way for heavy-light systems. This is done at the level of the 
interaction kernel eq. (|^) by simply allowing £ {u, d, s, c, 6} in (0). It has to be stressed that this extension is 
purely phenomenologicaly motivated, since the derivation leading to (0) is strictly valid for light quarks only. Also it 
should be mentioned that this use of the instanton induced interaction allows the mixing of heavy and light flavours 
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in the scalar and pseudoscalar flavour-neutral sectors. This has been investigated in || but will be neglected in this 
paper, since the effect turned out to be very small. 

Due to the J = selection rule of the residual interaction, a gradual fit of the model parameters to the light meson 
masses has been performed as explained in , leading to the parameter sets in table B. The resulting spectra of the 
isovector mesons are shown in fig. ([[]). The Regge trajectory is reproduced excellently in both models. The spectra of 
the light isoscalars and kaons are of similar quality. In particular the position of r\ and rj can be accounted for. The 
most significant difference between model A and B are the masses of the scalar mesons, esp. the a (980) in fig.(Q). 
For a detailed discussion of these spectra we refer to [Q and || . 

Based on this parameter set the masses of the heavy quarks have been adjusted to the J 7^ heavy mesons. Finally 
the additional couplings gf i f j have been adjusted to reproduce the heavy-light pseudoscalar masses. We want to 
stress that in this way all heavy mesons except for the pseudoscalar D, D Sl B and B s are described by fitting the 
heavy quark masses only. The resulting spectra are shown in figs. We find good agreement for the heavy-light 

mesons in both models with small advantages in model B due to a larger spin-orbit splitting. Also the gross structure 
of the heavy quarkonia spectra can be reproduced, but fine and hypcrfine splittings come out to small. This is due 
to the missing of a substantial spin-orbit and spin-spin interaction, as provided e.g. by the One Gluon Exchange 
potential. The radial excitations of the T, on the other hand, are excellently reproduced in model B up to the 5s 
state. Thus we think that we have reached a good overall agreement with the experimental heavy masses and therefore 
have gained a good estimate for the heavy meson amplitudes. To test these amplitudes further we investigate the 
semileptonic decays of heavy to light mesons. The relevant current matrix elements are calculated in lowest order 
after the prescription of Mandelstam |fl2[ . In our formalism this leads in general to 



(5) 



with the vertex function T p (p) := [S[ (-j + p)] 1 X P {p) \S F {— \- + p)] 1 ■ In the instantaneous approximation the 
vertex function in the rcstframe of the meson can be easily reconstructed from the Salpeter amplitude according to 

r p b)lp =( M,o) = -i J0fV(p,p>Mp>). (6) 

Thus the full 4-dimensional vertex function is known in the rest frame. Formal covariance of our model then allows 
to calculate the vertex function of any meson on its mass shell, which is essential for the calculation of meson decays. 



III. SEMILEPTONIC DECAYS 



The effective Lagrangian for the semileptonic decays, e.g. b — > c transitions, after integrating out the W boson, has 
the usual V-A current-current form 

C% = Vacrfil - 7 5 )& Ml " 1> (7) 

with the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V c &. Whereas the matrix elements of the leptonic 
current can be calculated exactly, those of the vector (V) and axial vector (A^) hadronic currents are parametrized 
by form factors, reflecting the hadronic structure. A common parametrization is Jl3| 

for a CP — ► 0~ transition, e.g. B — > Dlvi 

{D\V»\B) = f + (q 2 )(P£ + Pg) + f-{q 2 ){P" B - Pg) (8) 



and a — ► 1 transition, e.g. B — > D*£vi 



{D*\V\B) = J s>i ptT e*»prp^V( q 2 ) (9) 
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and 



(D*\A»\B) = 2m D « £ —^ q »A (q 2 ) + (m B +m D .) fe*" - ^Vj Mf) 

(n + ~ ^S^V) A 2 (g 2 ), (10) 

where q = P B — Pdm ' IS the momentum transfer and e the polarization vector of the vector meson. In the limit of 
vanishing lepton mass only 4 of these 6 form factors contribute to the decay rates, namely f+,V,A\ and A2, since 
terms proportional to <j M can be neglected. Then the differential decay rates, expressed in terms of these form factors, 
are 

T = '^ 2 ||l^(l^l 3 + ' ff -' 2 + l*' 2 >< < 12 » 

where we have introduced the hclicity amplitudes 

H ± (q 2 ) := (m B + mn^A^q 2 ) T 2 "^' V (q 2 ) (13) 
Ho{q2) ;= < -rnl.- q 2 imB+mD ^ Mq2) 2m|P*. A . 2 {q 2 ). (14) 



2m D ,^/q 2 m D *(m B + m D ,)^q 

With regard to the processes B — > D^iD, that is the transition between heavy flavours, we want to summarize the 
predictions of heavy-flavour symmetry. For these transitions, the Heavy Quark Effective Theorie (HQET) provides 
the appropriate framework. It is based on a systematic expansion in the inverse quark mass and has been worked 
out by Isgur and Wise jl4|,[l5). In the limit tuq — ► 00 the theory will become independent of the heavy degrees of 
freedom. For the spectra of heavy hadrons this will lead to degenerated doublets, corresponding to the two possible 
alignments of the heavy quark spin. Recent experimental results show evidence for this degeneracy in the spectrum 
of the D mesons For semileptonic decays HQET predicts a connection between the form factors of 0~ — > 0~ and 
0~ — > 1~ transitions. In particular, in the infinite quark mass limit, the transition matrix elements can be expressed 
by a single function only, the Isgur- Wise function £, which is normalized to unity at maximum recoil. 

In this heavy-quark limit it is more natural to express the decay amplitudes in terms of velocities rather than 
momenta and to introduce the dimensionless variable ui := UbUdoj = {m 2 B + m 2 DM — q 2 ) / (2mBm C (.)) instead of 
the momentum transfer q 2 . Therefore, a new set of form factors is used ]13J, defined by 



0- -> CT 



0- -> 1- 



{Dlv " lB) = h + {u){v B + «£) + M^)« - «£); (is) 



, = lZupa £ V B V D* h v {uj), 

(D*\A»\B) 



^m B m D , 



* M (w + 1) h Al (uj) - v B e* ■ v D * h A2 (uj) - v^,e* ■ v B h A3 (uj). (16) 



with 



uiq ^00 : h v = h Al = h As = h + = £ (17) 
h A2 = h- = 

in the infinite quark mass limit. The differential decay rates are given in this framework most conveniently by 
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dr 0_ 



4-0- 



do; 



^m%(m B +m D ) 2 (^ - lf"\V cb \^{u) 
G 



48tt 



^m 3 D »(m B - m D .) 2 y/uj 2 - l(w + if 



4u> 1 - 2wr* + r'l 
cj + 1 (l-r») 2 



|V cb | 2 ^ 2 ( W ) 



(18) 
(19) 



with r* = mo* /tub, where the two form factors Q and T are functions of h± and hy^h^i respectively. This 
parametrization has the advantage that for tuq — > oo J 7 and £7 become equal and coincide with the Isgur-Wise 
function £. 

In the following sections our results for semileptonic B and charmed meson decays are compared to the available 
experimental data and to the results of other models: the relativised constituent quark model of Isgur and Scora 
(ISGW2 fl6| ), the relativistic calculation of Wirbcl, Stech and Bauer in the infinite-momentum frame (WSB [ jl9| ) 
and the relativistic dispersion relation approach of Melikhov and Stech (MS [^0| ) . Since the experiments are not able 
to extract form factors from their measurements due to missing statistics, these are usually parametrized according 
to theoretical predictions. Therefore we find it convenient to parameterize our results in the same way, allowing a 
comparison with experimental data. The accuracy of these parametrizations is always indicated. 



A. Semileptonic decays of B mesons 



The decays B — > Divi and B — > D*£i/£ have been measured by Cleo and, more recently, by Opal |3|] 

and Cleo |§]. Figs, (g) and (§) show our results for the differential decay rate compared to the Cleo data 
p5|,53. We find good overall agreement with the experimental data for both decays, using a CKM matrix element of 
\V c b\ = 0.034 ± 0.001 and \V c b\ = 0.035 ± 0.001 for model A and B respectively, which has been determined by a x 2 
fit. These values are somewhat smaller than the PDG average of \V cb \ = 0.037 to 0.043. The resulting decay rates are 



model .A: T(B — ► DlVi) = 1.22 • 10 10 s -1 , T(B — > D*lvA = 3.21 • 10 10 " -1 



•s 



model B: T(B -> D£u e ) = 1.14 • 10 10 s -1 , T(B -> D*lv t ) = 3.24 • 10 10 -" 1 
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The comparison with the current world average of the Particle Data Group |lj 

PDG : T(B + -► D a t r v t ) = 1.30 ±0.13- 10 10 s _1 , T(B + -> D M i + v e ) = 3.21 ± 0.48 • lO^s" 1 
r(5° D'tvi) = 1.36 ± 0.12 • 10 10 s _1 , r(B° -> D*'t + vi) = 2.97 ± 0.17 • 10 10 s _1 



as well as the recent data of Cleo |2q] and Opal |32 



Cleo : T(B° -> D*+rv t ) = 3.66 ± 0.18 ± 0.23 • lO^s" 1 
Opal : T(B° -> D*+n? e ) = 3.40 ± 0.13 ± 0.30 • lO^s" 1 

shows satisfying agreement. But the new data tends to higher values of T, which would lead to higher values of \V c b\ 
from our calculation. 

In table [v| our results for polarization ratios, defined by 

T L _ fdq 2 q 2 P D *\H (q 2 )\ 2 T + _ J dq 2 q 2 P D , \H + (q 2 )\ 2 



T T fdq 2 q 2 P D *(\H + (q*)\* + \H_(q*)\zy T_ / dq*q 2 P D . \H_(q 2 



(20) 



are summarized. In these (| V c b\ independent) quantities we find good agreement with the experimental values as well. 
We also show the results for the form factor ratios 

= hy{ul = f q 2 \ V(ql) 

h A M \ (m B +m D *) 2 ) A 1 {q 2 ) 1 ' 

h Al (ui) V (n» B + m D *) 2 ) A x {q 2 ) 

In the heavy quark limit these ratios are predicted to be constant and equal to unity. In reality, due to corrections 
to this limit, R\ and i?2 do depend on lu, but at the scale of the 6-quark mass this dependency is expected to be 
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very weak. Therefore these ratios are prefcred in the analysis of B — > D*£v, where they are assumed to be constant, 
whereas h/n is approximated by a linear or quadratic function 

h Al H « h Al (1) (1 - p% (« - 1) + A Al (w - l) 2 ) . (23) 

For our form factors the w-dependency of R\/2 is less than 2%, and the mean values are shown in table M. The form 
factor h Al can be described by the parametrization (|2^) with an accuracy of 0.1% (0.3%) for model A(B), yielding 

model A: (1) = 0.97 ^=0.73 X Al = 0.27 

model B : h Al {\) = 1.01 p% = 0.98 A Al = 0.48, 

which is in good agreement with the experimental data shown in table [v|. 

Two recent measurements of B — > D*lv by the Cleo |28| and the Opal |32] collaborations have been published, 
where a different parametrization of the form factor based on dispersion relations has been used. These analyses 
expand hA 1 in the variable z{ui) = (\fu) + 1 — \/2~) /(•y/uT+T + \/2) and use the parametrization pl[ 

^(w)/^^) = 1 - 8p 2 z + (53p 2 - 15)z 2 - (231p 2 - 91)z 3 

with the (inconsistent) results 

Opal @ : p 2 = 1.21 ± 0.12 ± 0.20, Cleo @ : p 2 = 1.67 ± 0.11 ± 0.22. 

A similar fit to our form factor gives 

model A : p 2 = 0.83, model B : p 2 = 1.06 

with an accuracy of 1% and 0.4% for model A and B, respectively. Thus our calculation is rather compatible with 
the Opal result. 

To connect our results with the description in the framework of HQET, we have calculated the form factors T and 
G (see eq. (|l9|)) and their slopes at lo — 1 by fitting a quadratic function analogous to (|23|). This is possible with an 
accuracy of 0.2% (0.3%) for model A(B) and gives the result 

model A: Q{1) = 1.01, p% = 0.70, \ g = 0.22; ^(1) = 0.97, p\ = 0.65, Ajr = 0.20 
model 6 : 5(1) = 1.01, p% = 0.85, \g = 0.31; = 1.01, p 2 T = 0.91, A^- = 0.56 

These quantities 5(1) and T{V) have been calculated in the HQET. Current values are || 

5(1) = 1.00 ±0.07, T{\) = 0.92 ± 0.05. 

Finally we have performed the heavy quark limit numerically in our model by scaling the quark masses with a large 
factor. We find that the form factors then indeed coincide or vanish (see eq. (|l7|)). The resulting universal function 
£, which we identify as the Isgur-Wise function of our model, can be parametrized up to deviations of less than 0.2% 
and 0.5% for model A and B as 

£H = f(l)(l - p 2 (u; - 1) + A|(w - l) 2 ), 

where 

model A : f(l) = 1.00 pj = 0.78 A ? ~ = 0.31 
model B : f(l) = 1.00 p| = 1.06 A^ = 0.50. 

In particular, we find that £ is indeed normalised to £(1) = 1. The w-dependence is of course model dependent. 
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B. Semileptonic decays of charmed mesons 



The semileptonic decays of charmed mesons, induced by the quark level process c — > s, have been measured for the 



D — > as well as for the D s — ► rj/rj' /4> transitions. The results, averaged over isospin, are shown in table VI and 



VII| , For the D meson decay D — > Kiv the comparison with our calculation shows reasonable agreement for both 
models. The decay to the K* final state however is overestimated by about a factor of 2. The polarization obscrvables 
on the other hand are comparable with the experimental result, whereby model A gives better agreement than model 
B. This is a well known problem of constituent quark models. Our results are rather comparable with those of Wirbcl 



et al. (WSB |19|). With respect to the ISGW2 ||16[ predictions it is interesting to note that, whereas the inclusion 
of relativistic corrections was one of the main incredients in their model to decrease the D — > K* decay rate, this 
problem still exists in our full relativistic calculation. 

Apart from the decay rates and the polarization observables, the form factor ratios at zero momentum transfer, defined 
by 

n°) ^2(0) 
rv = mo)' r2 := awv (24) 

are considered. To extract these ratios from experiment the form factors are usually parametrized by a simple pole 
ansatz with a pole mass of 2.1 GeV for the vector form factor and 2.5 GeV for the axial form factors. We have 
performed such a fit to our calculations, which works up to deviations of about 6% for the form factors and of 4% for 
their ratios, and have extracted the form factor ratios for D meson decays from these parametrizations. The results 



in tabic VI show that the axial form factors arc generally overrated, while our vector form factor parameters agree 
with the values extracted from experiment. Thus the failure in the D — ► decay rate can be traced back to this 
overestimation of the axial form factors. 

In this connection it is interesting to estimate the effect of the full relativistic treatment. To do this we neglect 
the negative energy or 'lower' components of the Salpeter amplitude when reconstructing the vertex function. This 
procedure changes the norm of the Salpeter amplitude by a factor of 

m " °- 55 

for the D meson in model A. Here $++ := 14 2 7 (fr 1 ^ 7 and (.|.) is the scalar product induced by the Salpeter 
equation. This shows the importance of the negative energy components for heavy-light mesons. With this reduced 
vertex function the form factor ratios ry and r 2 are calculated. Hereby wo concentrate on model A, since, as has been 
discussed in || , the positive and negative energy components of the Salpeter amplitude decouple in the non-relativistic 
reduction of the Salpeter equation. We find that these ratios slightly rize: ry changes from 1.48 to 1.55, r^ goes from 
0.78 up to 0.840, while the q 2 dependence is hardly changed. Thus the modification of a form factors by the full 
relativistic treatment seems to be more complex than the intuitive correction anticipated in m& . 



The results on the semileptonic decays D s — > r\ and D s — > -q' arc shown in table VII. Here it has to be stressed that 
the flavour mixing of 77 and rj had already been fixed by the mass fit. No additional mixing parameter is necessary 
Although the differences between the results of our two models are rather large, the experimental data do not allow 
to prefer one of our parameter sets. 

Our results on D s — ► <f> show the same behaviour as those on D — > K*: While the decay rate is overestimated by a 
factor of 2, the polarization observable is in good agreement with the experimental data and the form factor ratio 
ry tends to be too small. Concerning the ratio of axial form factors it is interesting to note that the experimental 
value of r2 is 4 standard deviations higher than the corresponding value for D —> K* , which is in contradiction to an 
(approximate) flavour SU(3) symmetry. This result, if it should be confirmed, is clearly a challenge for any theoretical 
description. 

Finally, to conclude the discussion of charmed meson decays, we report our results on the Cabibbo suppressed process 



c — > d in tables VIII -X. The form factors again have been parametrized by a pole ansatz with pole masses of 2.0 GeV 
and 2.4 GeV for the vector and axial form factors respectively, which works up to deviations of 8% for Ai(0) and 5% 
for the ratios. The current experimental situation however allows no evaluation of our description of these decays. 



2 These are the calulated values, not to results of the pole ansatz parametrization in table |vi| 
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Thus in summary we find excellent agreement in the description of heavy to heavy transitions B — > £)M over the 
whole kinematic regime. We find also consistency with the description of these processes in the framework of the 
HQET. The results on heavy to light transitions are mostly in agreement with the experimental data, but the common 
problem of quark models to overestimate the axial form factors is also present here. 



IV. NON-LEPTONIC WEAK DECAYS 



To extract further information from our meson amplitudes we finally investigate non-leptonic decays. On tree level, 
non-leptonic decays arc mediated by a single W-boson emission. Hard and soft gluonic effects however might play a 
significant role in these processes. These corrections have been calculated with great effort in the last years in order to 
extract model independent Cabibbo-Kobayashi-Maskawa matrix elements or signatures for CP violation in B decays 
from experimental data. Thereby the hard and soft gluon contributions are separated by a Wilson operator product 
expansion, which results in the effective Lagrangian, e.g. for B — ► Dtt transitions, 

Zb^Ui = -^=V c bVj u Q (a 1 {fi){cb)^_ A {du)^ v _ A + a 2 (v){db)^_ A {cu)^ v _ A ^ 

+ l -C l ^){d\ A bY v _ A {c\ A u) ilv _ A + ±C 2 (p)(c\ A br v _ A (d\ A u), v _ A ^J (25) 

with (c\ A b)y_ A := cy(l — r y 5 )X A b etc., a\ := C\ + ^C 2 and a 2 := C 2 + 5C1, where C\/ 2 are the (scale dependent) 
Wilson coefficients, and X A the SU(3) color Gell-Mann matrices. Here we found it convenient to use an expression 
symmetric under Fierz transformation. The second line obviously does not contribute for color singlet states. 

This lagrangian gives rise to W emission diagrams. Contributions by weak annihilation and internal W exchange, 
which are suppressed by powers of Aqcd / n^b, are neglected. 

Thus the matrix element of a product of currents has to be evaluated. This is usually done using the "factorization 
approximation", where one assumes that the amplitude is dominated by its factorizable part. Then it is given by the 

product of two current matrix elements, e.g. for the transition B + — > D tt + 

A(B+^T?n+) = ^^^{ fll ( 7 r+|/ Vd jO)(^V^jB + >+«2(^Vv rfc |0>( 7 r + l^ b | J B + >}. (26) 

In this way the decay amplitude can be expressed by the decay constant and a form factor of the semileptonic decay 
at the relevant q 2 , e.g. 

(n+\J^\0)(D^\Jl\B + ) = (m 2 B m 2 D ) /.F (m 2 ) (27) 
<7r+|J^|0)<^°|j;|S + ) = 2e*-p B m D ,UA (ml) (28) 



with the decay constants (0| J^lO - ) = (0| J^l - ) = mFs^. It should be noted that the relevant form factors 

are -Fo and Aq, which arc unimportant for semileptonic decays due to the smallncss of the lepton mass. Hence non- 
lcptonic decays offer a possibility to access the remaining semileptonic form factors, at least within the framework of 
the factorization approximation. The factorization assumption has been proven recently at two loop order [ fl8f for a 
special class of decays. 

In this approach, however, the resulting amplitudes are scale dependent due to the /i-dependence of the Wilson 
coefficients, which is not canceled by the scale independent matrix elements. To deal with this difficulty the coefficients 
ai and a 2 are often treated as effective and free parameters, corresponding to some factorization scale, to be extracted 
from the data. But since we are interested in an estimate of the quality of our form factors, we find it sufficient to 
neglect all strong gluonic effects, which results in C\ = 1, C 2 = = a\ = 1, a 2 = 5 and restrict our calculation to the 
decays of B mesons via the emission of a W meson, usually called Type I decays. 



The weak decay constants are shown in table Kill. Since these are generally overestimated in our models we have 



used the experimental values, which are extracted from leptonic decay or t— decay and are summarized in table XIII 



For the vector mesons D* and D*, where no data are available yet, we use Fd* ~ fD (a) , which is valid in the heavy 



quark limit. Deviations from this limit are expected to be about 10-20 % [|22|. The relevant CKM-matrix elements 
are taken from H except for V c b where we use the fit results from section III A . 
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Our results for non-leptonic B-dccays are shown in table XII , compared with the experimental data from M as well 
as the calculation of Neubert et al. []22f] . We find good agreement with the data available so far for both our models; 
and for those decays, which are not measured yet, our results are comparable with pi| . 

In order to stress the influence of our form factors we show the ratios of decay rates (in which the decay constants as 
well as the coefficient a\ cancel) in table XI. We also find good agreement with the experimental data, however the 
errors are quite large. 



V. SUMMARY 



In this paper we have calculated the masses and the exclusive semileptonic and non-leptonic decays of heavy mesons in 
a constituent quark model based on the Bcthe-Salpeter equation in instantaneous approximation. Apart from a linear 
confinement potential with two different phenomenological Dirac structures, a flavour dependent residual interaction 
motivated by instanton effects is adopted, which has been naively generalized to heavy flavours. Thus extending a 
very good description of light mesons, which has recently been updated in [Q, we also find good overall agreement 
with the data on heavy meson masses. 

The resulting meson amplitudes are used to calculate semileptonic decay rates. We find excellent agreement in the 
description of heavy to heavy transitions B — > over the whole kinematic regime. Our results are also consistent 
with the description of these processes in the framework of the HQET. The results on heavy to light transitions are 
mostly in agreement with the experimental data, but the common problem of quark models to overestimate the axial 
form factors is also found here. 

Finally non-leptonic decays have then been calculated in the approximation of factorizing matrix elements. In spite 
of this simple picture we find good agreement with the experimental results on B meson decays. 

Acknowledgement: Financial support by funds provided by the Graduicrtcnkolleg "Die Erforschung subnuklearer 
Strukturen der Materie" is gratefully acknowledged. 
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FIG. 1. The spectrum of the isovector mesons. The centric column shows the experimental data from Errors are indicated 
by shadowed boxes. The left and right columns show our results with model A and B, respectively. 
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FIG. 2. The spectra of the charmed D and D a mesons. The centric column shows the experimental data from [1]. Errors 
are indicated by shadowed boxes. The levels marked with t are taken from [Q. The left and right columns show our results 
with model A and B, respectively. Note that the 1 + states are twofold degenerated in our calculation corresponding to the 
total spin S — and S = 1. 
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FIG. 3. The spectra of the bottom B and B a mesons. The centric column shows the experimental data from [[fj. Errors are 
indicated by shadowed boxes. The levels marked with t are taken from [Q. The left and right columns show our results with 
model A and B, respectively. Note that the 1 + states are twofold degenerated in our calculation corresponding to the total 
spin S — and S = 1. 
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FIG. 4. The spectra of heavy quarkonia. The centric column shows the experimental data from 
shadowed boxes. The left and right columns show our results with model A and £>, respectively. 
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FIG. 5. The differential decay rate for B — > D£u. The data are taken from The values of \V c b\ corresponding to our 

calculations are |V c i,| = 0.034 and |V c b| = 0.035 for model .4 and £>, respectively. 
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TABLES 



Parameter Model A Model B 





m„ 


[MeV] 


306 


380 


Constituent 


m s 


[MeV] 


503 


550 


quark masses 


m c 


[MeV] 


1835 


1780 




i>H, 


[MeV] 


5240 


5092 




Qnn 


[GeV" 2 ] 


1.73 


1.62 




9 ns 


[GeV- 2 ] 


1.54 


1.35 


Residual 


Que 


[GeV" 2 ] 


1.11 


1.58 


interaction 


9nb 


[GeV" 2 ] 


0.53 


1.07 




g sc 


[GeV" 2 ] 


0.65 


1.27 




9sb 


[GeV" 2 ] 


0.00 


0.76 




A 


H 


0.30 


0.42 


Confinement 


a c 


[GeV] 


-1.751 


-1.135 


parameters 


b c 


[GeV fm- 1 ] 


2.076 


1.300 


Spin structure 




or 


\ (l ® 1 -7° ® 7 °) 


1 (l ® 1 - 7 s <g> 7 5 - 



7m ® 7 



TABLE I. The parameters of the confinement force, the 't Hooft interaction and the constituent quark masses in the models 
A and B. 
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Meson {.F) 


n 


Model A 


Model B 


Meson (J") 


n 


Model A 


Model B 


7-1 / n — \ 

D(0 ) 





1869 


1869 


7~i fr>— \ 
L>s (0 ) 





1969 


1969 




1 


2677 


2578 




1 


2794 


2683 




2 


3258 


3041 




2 


3388 


3152 


zr(i-) 





1993 


2034 


D* S {1~) 





2049 


2116 




1 


2769 


2648 




1 


2852 


2746 




2 


2822 


2679 




2 


2905 


2776 




3 


3327 


3079 




3 


3432 


3192 




4 


3344 


3101 




4 


A r A 

3454 


3210 


7-1* / r\-\- \ 





2519 


2375 


D a0 {0 ) 





2563 


2464 




1 


3115 


2884 




1 


3196 


2986 


D i{l ) 





2464 


2420 







2532 


2506 




1 


2464 


2420 




1 


2532 


2506 




2 


3075 


2908 




2 


3172 


3007 




3 


3075 


2908 




3 


3172 


3007 


D* 2 {2+) 





2475 


2469 


D* s2 (2 + ) 





2541 


2552 




1 


3086 


2930 




1 


3182 


3032 




2 


3128 


2957 




2 


3223 


3057 


TABLE II. 


Masses of the charmed D and D 3 


mesons in MeV, calculated in model A and B; 


n denotes the radial excitation. 


Note that the 1 + states are twofold degenerated correspondin 


g to the total spin S 


= and S -- 


= 1. 




Meson (J n ) 


n 


Model A 


Model B 


Meson (J n ) 


n 


Model A 


Model B 


B(O-) 





5279 


5279 


B s (0~) 





5368 


5369 




1 


6002 


5869 




1 


6101 


5960 




2 


6517 


6279 




2 


6628 


6376 


B*(l~) 





5325 


5346 


b;(i~) 





5369 


5425 




1 


6035 


5900 




1 


6102 


5986 




2 


6108 


5947 




2 


6172 


6031 




3 


6543 


6296 




3 


6629 


6392 




4 


6575 


6322 




4 


6664 


6417 


Bo(0+) 





5796 


5675 


B* a0 (0 + ) 





5822 


5750 




1 


6341 


6120 




1 


6401 


6208 


Bi(l + ) 





5770 


5696 


B sl (l+) 





5822 


5774 




1 


5770 


5696 




1 


5822 


5774 




2 


6322 


6136 




2 


6401 


6224 




3 


6322 


6136 




3 


6401 


6224 


B 2 *(2+) 





5771 


5711 


BU2+) 





5823 


5788 




1 


6324 


6140 




1 


6402 


6230 




2 


6391 


6186 




2 


6466 


6273 



TABLE III. Masses of the bottom B and B s mesons in MeV, calculated in model A and B; n denotes the radial excitation. 
Note that the 1 + states are twofold degenerated corresponding to the total spin 5* = and 5 = 1. 
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IVltJoOIl 1 J J 


11 


Model A Model B 


ivieson i j i ii 




iviouei o 







3037 


3114 




Vb(0~ + ) o 


9497 


9562 




1 


3789 


3708 




1 


10055 


9983 




2 


4363 


4147 




2 


10503 


10318 


J/*(l— ) 





3039 


9 1 9 9 

3133 




T(l~ ) 


9497 


9565 




1 


3792 


3719 




1 


10055 


9985 




2 


3846 


3754 




2 


10103 


10020 




O 
O 


4366 


4155 




O 
O 


iuouo 






4 


4398 


4176 




4 


10534 


10342 












5 


10892 


10609 

















1UDZO 












7 


11250 


10871 












8 


11261 


10882 


h c (l+-) 





3482 


3483 




Mi+-) 


9824 


9811 




1 


4110 


3965 




1 


10304 


10171 


Xco(0 ++ ) 





3482 


3449 




Xto(0 ++ ) 


9824 


9805 




1 


4110 


3943 




1 


10304 


10166 







3482 


3483 




Xi.i(l ++ ) o 


9824 


9811 




1 


4110 


3965 




1 


10304 


10171 


Xc2(2++) 





3485 


3508 




Xb2(2++) 


9824 


9816 




1 


4114 


3983 




1 


10304 


10175 




2 


4164 


4015 




2 


10352 


10210 


TABLE IV. Masses of the heavy quarkonia 


in MeV, calculated in model A and B; n denotes the radial excitation. 


Parameter 




Exp. § 


Model A Model B ISGW2 [[[§ 


WSB 


MS @ 


T(B — » D) [l0 13 |V ci) 








1.05 


0.93 1.19 


0.808 


0.86 


r(B — » D*) [lO^IVcbl'a -1 ] 






2.78 


2.64 2.48 


2.19 


2.28 


r L /r T 




1.24 ±0.16 @ 




1.14 


1.20 1.04 




1.11 


r+/r_ 








0.23 


0.27 






Ri 




1.18 ±0.30 ±0.12 




1.18 


1.10 1.27 


1.09* 




R> 




0.71 ±0.22 ±0.07 




0.94 


0.87 1.02 


1.06* 




2 

Pa 1 




0.91 ±0.15 ±0.06 




0.75 


1.02 







TABLE V. B — > D^l v decay observables and form factor parameters. We use \V c i 
described in section III A. 



0.034(.4) resp. |V c i,| 



* taken from jlj 

0.035(B), as 
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Parameter 



Exp. §] Model A Model B ISGW2 |l6| 



WSB ^ MS |§ 



r(D — > A) 1 10 s \ 


7.97 ± 0.36 


7.51 


7.26 


10.0 


8.26 


9.7 


F(D — > A*) [io^s- 1 ] 


4.55 ±0.34 


7.64 


10.08 


5.4 


9.53 


6.0 


Tl/Tt 


1.14 ±0.08 


1.29 


1.48 


0.94 


0.91 


1.28 


r+/r_ 


0.21 ±0.04 


0.23 


0.34 








Ai (01 


0.56 ±0.04 H 


0.69 


0.81 




0.88 


0.66 


ry 


1.82 ±0.09 


1.54 


1.18 


2.0* 


1.44 


1.56 


''2 


0.78 ±0.07 


0.81 


0.62 


1.3* 


1.31 


0.74 


TABLE VI. D -> A' 
isospin. We use |14 s | = 


decay observables and form factor parameters. 
0.975 §. 


The experimental decay rates are 


* taken from j^sj] 

averaged over 


Parameter 


Exp. § 


Model A 


Model B 


ISGW2 [16 


WSB @ 


MS @ 


T(D S ^ V ) [lO 10 ^ 1 ] 


5.24 ± 1.41 


4.05 


3.11 


3.5 a) 




5.0 


V(D 3 -> 77') [I0 lu s L J 


1.80 ±0.69 


1.27 


1.75 


3.0 a) 




1.85 


r(£>.->^) [io 10 ^ 1 ] 


4.03 ± 1.01 


7.89 


9.67 


4.6 




5.1 




0.72 ±0.18 


1.20 


1.42 


0.96 






r+/r_ 




0.20 


0.33 








Ai(0) 




0.66 


0.79 






0.65 


rv 


1.92 ±0.32 


1.77 


1.30 


2.1* 




1.71 


r2 


1.60 ±0.24 


0.85 


0.63 


1.3* 




0.72 







TABLE VII. D 3 — > rj/rj' /(f>£u decay observables and form factor parameters. The experimental decay rates are averaged over 
isospin. We use |V CS | = 0.975 §. 

A 77 — rf mixing angle of —20° is assumed. An angle of —10° would lead to 5.3 and 2.3. 
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Parameter 


Exp. § 


Model .4 


Model B 


ISGW2 @ 


WSB JlJ 


| MS @ 


T(D + ^7v°) [lO 9 *- 1 ] 


2.9 ± 1.4 


1.03 


0.99 


2.4 


3.6 


4.8 


r(D°^n-) [lO 9 ^ 1 ] 


9.0 ±1.5 


2.06 


1.99 


4.8 


7.1 


9.6 


T (D + ^p°) [io 9 ^ 1 ] 


2.1 ±0.8 


2.25 


3.36 


1.2 


3.4 


2.1 


Tl/Tt 




1.30 


1.55 


0.67 


0.91 


1.16 


r+/r_ 




0.15 


0.26 










— 


0.58 


0.72 




0.78 


0.6 


ry 




1.68 


1.26 




1.58 


1.48 


''2 


— 


0.77 


0.59 




1.18 


0.82 


TABLE VIII. 


D -> Tr/pte 


decay observables and form factor parameters. We use 


|V cd | =0.222 




Parameter 


Exp. § 


Model A 


Model B 


ISGW2 


WSB [19 


MS @ 




— 


0.79 


0.95 


1.5 a) 


— 


— 


r(D^r ? ') [lo'^- 1 ] 


— 


0.14 


0.16 


0.3 a) 


— 


— 


r(£>->w) [io^- 1 ] 




2.23 


3.35 


1.2 






Tl/Tt 


— 


1.31 


1.55 


0.68 






T+/Y- 




0.15 


0.26 








Ai(0) 




0.41 


0.51 








ry 




1.68 


1.27 








r2 




0.77 


0.59 









TABLE IX. D — > 77/77' /totv decay observables and form factor parameters. We use \V c d\ = 0.222 |l| 
A 77 — 77' mixing angle of —20° is assumed. An angle of —10° would lead to 1.1 and 0.4. 



Parameter 


Exp. § 


Model .4 


Model B 


ISGW2 @ 


WSB JT§ 


MS @ 


Y{D S ^K) [lO 9 ^ 1 ] 




3.42 


3.15 


4.4 




6.4 


Y{D a ^K") [lO 9 *- 1 ] 




2.71 


4.54 


2.2 




3.9 


Tl/Tt 




1.24 


1.51 


0.76 




1.21 


r+/r_ 




0.14 


0.29 








Ai(0) 




0.43 


0.58 






0.57 


ry 




1.86 


1.26 






1.83 


r2 




0.76 


0.56 






0.74 



TABLE X. D s —* K/K*lv decay observables and form factor parameters. We use \V c d\ = 0.222 Jl). 
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Decay ratio 



Exp. 



Model A 



Model B 



NRSX 



r(B°-.D~ir+) 



r(B u ->i5* 


-7T+) 


r(s°— >D _ 


-P + ) 


r(s°->£)* 


-P+) 


r(s°— >D" 




r(s°^i5* 




r(B+- > S°£)+) 


r(B+^D* 


V) 


r(s°— ►£)" 




r(s°^i5* 









T(S+^D £>* + ) 



D~Df 

D~D* S + 
>D*°Df 



B+ - 
B° -» 
B+ - 
5° -» 
B+ - 

b+ -» 

b () d*~d: 



1.09 ±0.17 
1.16 ±0.62 

0.83 ±0.43 
1.08 ±0.56 

0.50 ±0.31 
0.33 ±0.19 



1.00 
0.87 

1.42 
0.41 



1.05 
0.92 

1.28 
0.36 



7.9 ±2.4 

5.2 ± 1.9 

5.4 ±2.4 

6.5 ±3.2 

7.3 ±3.0 
6.2 ± 2.2 
16.3 ±6.0 
12.9 ±4.5 



8.73 
8.73 
6.05 
6.05 
6.14 
6.14 
14.68 
14.68 



8.27 
8.27 
5.49 
5.49 
6.47 
6.47 
15.24 
15.24 



1.04 
0.88 

1.47 
0.39 





TABLE XL 


Ratios of non-leptonic B meson decay rates 




Decay mode 


Exp. § 


Model A Model B 


NRSX 



B° -» D~n+ 
B° -> D-p+ 
B° -^D*-7r+ 
B° -^D*-p+ 


1.94 ±0.26 
5.10 ±0.90 
1.78 ±0.14 
4.4 ±2.2 


2.21 
5.68 
2.22 
6.57 


1.97 
6.13 
1.88 
6.67 


1.94 
4.84 
1.87 
5.48 


B° -> D~K+ 
B° -> D~K* + 
B° -> D*-K+ 
B° -> 




0.17 
0.31 
0.16 
0.36 


0.15 
0.30 
0.14 
0.34 


0.15 
0.26 
0.14 
0.32 


5° -» 

5° -> D*~D + 
B° — ^ D*-D*+ 


-0.20 


0.53 
0.37 
0.39 
0.86 


0.49 
0.34 
0.41 
0.89 


0.23 
0.23 
0.17 
0.54 



6.61 
6.61 
6.15 
6.15 
4.49 
4.49 
15.8 
15.8 



TABLE XII. Non-leptonic B decay rates Y [n£_lL We use |V CS | = 0.975, \V ud \ = 0.975, \V ua \ = 0.223 § and \V cb \ = 0.034 (A) 
resp. \V c b\ = 0.035 (B), as described in section III A. The results of p3| have been adapted to the decay constants used in our 
calculation. 
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Meson 



Exp. § 



Model A 



Model B 



n 
K 
D 
D, 



130.7 ±0.1 ±0.36 

159.8 ± 1.4 ±0.44 

o nn +180+80 
OUU -150-40 

280 ± 19 ± 28 ± 34 



212 
248 
293 
315 



219 
238 
263 
284 



P 
u> 

K* 
D* 

d: 



216 ±5 1 22] 
195 ±3 [~] 

237 ±4 [22] 
214± 7 22] 



470 
472 
475 
467 
339 
378 



717 
726 
685 
695 
409 
445 



TABLE XIII. Decay constants of pseudoscalar and pseudovector mesons in MeV 
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